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Abstract 



We present a closed form expression for |V^&| 2 /|V#,V^*| 2 in terms of the 
endpoint photon and lepton spectra from the inclusive decays B — > X s 7 and 
B — > X u £v, respectively, which includes the resummation of the endpoint 
logs at next to leading order and is completely independent of the B meson 
structure function. The use of this expression for extracting V u b would elimi- 
nate the large systematic errors usually incurred due to the modeling of the 
heavy quarks' Fermi motion. 
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I. INTRODUCTION 



Presently great effort is being put into determining the nature of CP violation in the 
Standard Model. The unitarity triangle is being constrained with increasing precision via a 
combination of measurements. The size of the allowed regions in the (p, rf) plane |1| could be 
greatly reduced if we could more precisely calculate the hadronic dynamics involved in weak 
decays. Indeed, one of the largest uncertainties in determining the position of the upper 
vertex of the unitarity triangle stems from our poor knowledge of the CKM element V u b- 

Extracting V uh presents quite a challenge to the theoretical community. The elegant 
methods of heavy quark effective theory (HQET) used in determining V c b from exclusive 
B to D transitions are of little use in the heavy to light modes ||. Whereas the extraction 
from inclusive decays is hindered by the experimental cuts which force theorists to make pre- 
dictions in a region of phase space which is especially sensitive to non-perturbative hadronic 
dynamics. The presence of these experimental cuts thus introduces theoretical errors, due 
to the use of models, which are difficult, if not impossible, to quantify. Here we will present 
a method for extracting \V u b\ which has the benefit of avoiding the issue of the unknown 
hadronic dynamics, thus obviating the need for any models. 

To extract \V u b\ from inclusive semi-leptonic B decays, we may study the endpoint region 
of the charged lepton energy spectrum in B — > X u C,v decays and apply a lower cut on Ef 
above the b — > c endpoint energy in order to eliminate the background from charmed decays. 
In the Standard Model, this decay rate can be calculated in a systematic expansion in a s , 
1/ (Mw, mt) and l/m&. Large amounts of resources have gone into calculating the decay rate 
at next to leading order in the strong coupling. The calculation is broken into three stages. 
The full standard model is first matched onto the four fermion theory at the scale Mw- It 
is then run down to the scale m&, after which the inclusive decay rate may be calculated in 
HQET |4],[| in expansion in l/m&. 

However, near the endpoint region certain effects arise which jeopardize the accuracy of 
the calculation. The expansion in a s breaks down in the endpoint region due to presence 
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of the threshold logs, log(l — 2E^/m&). In addition, the expansion in l/m& breaks down 
because the relevant expansion parameter in this region is l/(m& — 2Eg) = l/[m&(l — x)] 
and the true kinematic endpoint is determined by the hadronic mass Mb, instead of the 
partonic mass m&. Thus, both perturbative and non-perturbative expansions break down at 
the endpoint. Given the proximity of the b — > u and b — > c endpoints, there is only a very 
narrow window, AE ~ 350 MeV, left for the extraction of \V u b\. In particular, due to the fact 
that there is an imposed energy cut on the charged lepton, the small parameter 1 — 2Ei/m\ ) 
entering the calculation hinders our ability to make reliable theoretical predictions. Both the 
perturbative expansion in a s as well as the non-perturbative expansion in must 
be reorganized, in a systematic expansion in powers of a s log(l — x) and AQ CD /[m b (l — x)], 
respectively. While the perturbative expansion is, in principle, calculable, the reorganiza- 
tion of the non-perturbative expansion necessitates the introduction of a non-perturbative 
structure function, of which we are ignorant. Thus, it would seem that we are left in a lurch 
as far as an extraction of \V u b\ from inclusive decays is concerned. 

However, as is well known |^J7|], the aforementioned structure function is universal in B 
decays, in that it is independent of the short distance decay under consideration. Therefore, 
we may extract it from another decay, such as B — > X^ for use in the prediction for 
B — > X u lv. However, given the dearth of data near the endpoint this extraction will incur 
large errors, given that we don't know its shape. It would be preferable if we could sidestep 
the extraction completely. In this paper we will show that the process of extraction of the 
structure function may indeed be bypassed, and that a prediction for the ratio |V^| 2 /| V^*V^| 2 
may be given at next to leading order in log(l — x) without model dependence. 



II. THE ONE LOOP RESULT 

The calculation of inclusive semi-leptonic B — > X u iv decays begins with the low-energy 
effective Hamiltonian |8| 

H eff = -^v ub (u lfl p L b) (I^Pm) 
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where Gp is the Fermi constant and is the left-handed projection operator |(1 — 75). 
J M and Jg are the hadronic and leptonic currents, respectively. The differential decay 
distribution is 

d 3 T 



dq 2 dE^dE v 

where q 2 = (pe + p u ) 2 and the leptonic tensor is 



2G 2 F \V vb \ 2 W aP L af} , 



L at3 = 2 (p?p£ + $pl - g a % ■ p v - i^ Xa p^p vX ). 



(2) 



(3) 



The hadronic tensor is related to the imaginary part of a time ordered product of hadronic 
currents 
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d A x e-^ x (B\T(Jl(x)MO))\B). 



(4) 
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It is this time ordered product which can be computed using an operator product expan- 
sion (OPE) in terms of local operators, involving 6-quark fields in HQET, whose Wilson 
coefficients can be calculated in perturbation theory away from the cuts of T Q/ g [|],[5|. In 
particular, the one loop decay spectrum, including the leading non-perturbative corrections, 
is given by M 
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G(x) = log 2 (l -x) + 2Li 2 (x) + — + 
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and x = 2Ei/mb. In Eq. (Q) the first line is the non-perturbative corrections to order 
(AgcD/^b) 2 and the second line contains the perturbative corrections to order a s . The 
non-perturbative corrections to this order are characterized by only two hadronic matrix 
elements 

\ 1 = {B(y)\b v {iDfb v \B(y)), (9) 
A 2 = \(B{v)\b v 9 -a^ u b v \B(v)). (10) 

Here b v is the velocity dependent 6-quark field as defined in HQET, g is the strong coupling, 
and is the strong interaction field strength tensor. A2 can be determined from the 
B — B* mass difference, yielding A2 ~ 0.12 GeV 2 . From the above expression we see that the 
breakdown of the expansions, in a s and AqcD/^b, manifest themselves in the large logs and 
the derivative of delta functions, respectively. Though the differential decay rate is singular 
as x — > 1, these endpoint singularities are integrable and the total decay rate is given by 
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The maximum lepton energy for a particular hadronic final state X is Ei = (M| — 
Mx)/(2Mb) in the B rest frame, and consequently leptons with energies in the endpoint 
region (neglecting the pion mass) Mb/2 > E% > (M B — M^) j '(2Mb) can arise only from 
b — ► u transition. This translates into a lepton energy cut x c ~ 0.87. 



III. THE PERTURBATIVE RESUMMATION 



The systematics of the perturbative expansion is discussed at length in Ref. fT0"| . Here 
we summarize the discussion. The perturbative series near the endpoint schematically takes 
the form 

Coo 

+ C 12 a s log 2 (l -x) + C u a s log(l - x) + C 10 a s 

+ C 2A a 2 s log 4 (l -x) + C 23 a 2 log 3 (l - x) + C 22 a 2 s log 2 (l - x) + C 21 a 2 s log(l - x) + C 20 a 2 s 



+ C 36 a 3 s log 6 (l -x) + C 35 ^log 5 (l -x) + C M a 3 \og\l -x) + C 33 a 3 \og 3 (l - x) + ■ 



+ 
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Naively, if we would like to go to a region where a™ log n (l — x) ~ 1, it seems that we would 
have to sum a large number of subleading logarithms below the line C m k, n k, a daunting task. 

However, this criteria is not correct. The reason for this is that the series is known to 
exponentiate into a particular form. In fact, the series resums into a function of the form 

~dT~ 



log 



dx 



log(l - x)g 1 [a 8 f3 log(l - x)\ + g 2 [a s (3 log(l - x)] 
+a s g 3 [a s p Q \og(l — x)] H , 



(13) 



where /3q = (33 — 2Nf) / (12tv). This result implies that the aforementioned triangle assumes 
a definite structure. The resummation does not sum the entire triangle, but the terms 
dropped are subleading in the exponent. For more details see [ [1C[ . 

The resummation itself is performed in moment space, where the rate factorizes into 
a short distance hard part (H N ), a soft part (ojv), and a jet function (J N ) |H| . Using 



renormalization group techniques it is possible to show that [12 



a N J N = exp [log Ngx(x) + ^(x)] 



(14) 



where x — a s( m b)Po\°gN, 



and 
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In the above, ft = (153 - 19 iV / )/(24vr 2 ), k = 67/6 - tt 2 /2 - 5JV//9, x v = 2E„/m b , and 
je — 0.577216... is the Euler-Mascheroni constant. The moment space resummation can be 



written as [|12|] 

1 f MB/mb ,^N-i d dT 



M N = -— dxx 1 

1 o -Jo ax ax 

= S N £ dx v 6(1 - x v )x v H{x v )e l ^ N9 ^ +9 ^ + O C^j . (18) 

Note that in writing the above we have taken into account moments of the non-perturbative 
structure function, Sjv, which arises from resumming leading twist terms in non-perturbative 
expansion in l/m& (see [|1^] for details.) The hard part is given byQ 



H(x u ) = l 



3tt 



2tt 2 5.0/-, 3 .. , . log(l — x v ) nT . , 
+ - + log 2 (l -x v )-- log(l - x„ + ^ + 2Li 2 (x, 



• (19) 



3 2 ° x 2 ov ' y 

To go back to x-space, we must take the inverse Mellin transform of Eq. ([Tj|) which is 
given by 

1 dT 1 fC+ioc dN 

r P~~f = 7T- / M N+1 X- N — 

1 o dx 2m Jc-ioo N 

rl l rC+ioo dN 

= dx u 6(l - x v )x v H{x v )— I —x~ N S N+1 a N+1 J N+ i. (20) 

JO 2,1X1 JC-ioo iV 

Then using the identity derived in the Appendix, and ignoring Sn, we obtain 

1 r C + io ° dN 



^x~ N exp{log(iV + 1) gi [a s (3 \og(N + 1)] + g 2 [a s (3 log(JV + 1)]} 

2711 JC-ioo iV 

x [1 + F(a„l)], (21) 



e l g 1 (a s /3 l)+g2{a s l3ol) 



r [1 - gi(a s f3 l) - otsPol g'x(a s /3 l)] 
where 



%,i) = W E/yP ( 22 ) 

k=l j=0 



x The hard part given here is different from that in Ref. [j!2[| . In the Appendix we perform the 
inverse Mellin transform to next-to-leading log accuracy, therefore the hard part can differ since it 
is a next-to-next-to-leading contribution. The hard part here is chosen to reproduce the one loop 
rate when x — > 1 after performing the inverse Mellin transform given in the Appendix. 
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represents next-to-next-to-leading log contributions and 

/ = — log (— logo;) 
~ — log(l — x) when x ~ 1. (23) 

It is important at this point to note a crucial difference between this calculation and 
resummations carried out in threshold production, which has been observed in B — > X s 7 
decays as well. In threshold production (Drell-Yan for instance), the overall sign of g\ is 
flipped (in standard schemes). This difference in sign changes the nature of the inverse Mellin 
transform. In the case of a negative exponent (our case), the inverse Mellin transform is 
integrable, when expanding g\ and in the Sudakov exponent, whereas it is not the case 
for threshold production. The expansion of the series for the Drell-Yan case at leading 
logarithms accuracy has introduced factorial growth terms at large order, which shows up 
as a non-integrable pole on the positive axis of the Borel plane. If we were to perform 
the inverse Mellin transform exactly, this spurious factorial growth would be absent, as was 
pointed out by Catani et al. |13|]. It would be interesting to apply the analytic formulae 
derived in the Appendix to study the behavior of these factorially growing terms. 

We may use these results to pick out certain terms in the two loop calculation. In 
particular, we may expand the analytic expression and determine the coefficients of the 
terms of order a;?log 4 (i — x) and a;?log 3 (l — x). Lower order logs will not be reproduced 
correctly given the fact that we have dropped terms of order a s in the exponent. Expanding 
we find 



r dx 



= ^ g \l-x) + ^(\^ + f_)\o g \l- ,■)-- (21.) 
e>(«2) J7r it \6 zt 



The term proportional to a^/3o is P ar t of the BLM correction, which is typically a good 
approximation to the full two loop result. However, note that the contribution of the last 
term in Eq. (pl|) , not proportional to (3q, is actually larger than the corresponding BLM 
correction, the second term in Eq. (£23|). Similar situation has been observed in B — > X s 7 



decays in Ref. [10]. Although one can not rule out the possibility of a cancellation with other 



terms such that the BLM result in the two loop result still dominates, this does cast doubt 
on the issue of BLM dominance in the endpoint region, since we know that Eq. ( p^) is the 
leading contribution at 0(a 2 s ) when x — > 1. Terms we neglect in Eq. (|24| ) are suppressed by 
at least a factor of log(l — x) in the endpoint region. 



IV. B -► X s -f DECAYS 

Here we briefly overview the inclusive B — ► decays and give the analytic ex- 
pression for the resummed rate. For B — ► X s 7 decays, only the high energy part of the 
photon spectrum is experimentally accessible due to large background cuts. Currently the 



experimental cut on the photon energy is given by E° ut > 2.1 GeV | 15fl . This introduces a 



small dimensionless parameter 1 — 2£^ Mi / m ft into the calculation, resulting in large threshold 
logs, log(l — 2E^ ut /mi > ), which grow parametrically as E^ ut approaches the endpoint. 

The calculation of inclusive decay rates begins with the low-energy effective Hamiltonian 

AC 8 

Heff = CiWOiiy), (25) 

V 2 i= i 

where Cj(//) are Wilson coefficients evaluated at a subtraction point /i, and Oj(/i) are di- 
mension six operators. For B — > X s 7, the only operators that give a relevant contribution 
are 

2 = (c La 'fb La ){s L p'y ll ,C L p), 
g 

°7 = -T^-^rribSLaCr^F^bna, (26) 



16tt 2 



°8 = TTr^mbSLaV^G^T^bnp. 



Near the endpoint, the decay rate is dominated by the Oj operator. In particular, the 
decay rate due to the Of operator is 

1 dT" 1 

*<i 71 3 (! " x ) 



Tq dx 
where x = 2E 1 /m\ ) and 



a s (2x 2 -3x- 6)z + 2{x 2 - 3) log(l - x) 

5 (27) 



As x — » 1 this contribution diverges as log(l — — x). The integrated rate is finite due 

to virtual corrections at x — 1. The differential rate in moment space factorizes into the 



form |[T 



a N Jjj = exp [log Ng x {x) + 9lix)\ , (29) 
where giix) is the same as in the semi-leptonic decays Eq. (|15D and g\ (x) = ^(x) - gstiXi x v) 



T2fl . #2(x) an( i gsiiXi x u) are defined in the semi-leptonic decays Eq. (|16|) and we see that 
the only difference between these two cases is that in the semi-leptonic decays we have an 
additional g s i piece in the Sudakov exponent and the extra x u integration. The moment 



space resummation is written as 12 



1 [ MB i mb , 

Ml = -= / dxx N ~ l 



v " Jo dx 

S N H~<e loeN9l{x)+g 2 {x \ (30) 



and the hard part is given by 



It is important to note here that the non-perturbative structure function Sn is universal in 
B decays and therefore the same as in the inclusive semi-leptonic decays. Evaluating the 
inverse Mellin transform using Eq. (|A12 ) and ignoring Sn, we get 



1 dT"< 1 /-C+ioo 

' ' 1 Mlx~ N dN 



Tq dx 2?ri Jc 



-too 

,1/ 



d C e l9i{a s /3 l)+g2(a s /3ol) ~| 

~ % Tx X (1 ~~ X)H '^ [1 - g~MM - a s (3 lg[(a s {3 l)} j ' (32) 



V. THE NON-PERTURBATIVE RESUMMATION 



Before we can confront theoretical calculations presented so far with experimental data, 
we have to take into account non-perturbative corrections. As is well known, as one 

10 



probes the spectrum closer the end point, the OPE breaks down, and the leading twist 
non-perturbative corrections must be resummed into the B meson structure function ||[7||. 
Formally we may write the light cone distribution function for the heavy quark inside the 
meson as 

f(k+) = (B(v) | b v 5(k+ - iD+)b v | B(v)). (33) 

While the shape of this function is unknown, the first few moments of f{k + ), 

A n = J dk + k n + f(k + ) 

= (B(v) | b v (iD + ) n b v | B(v)), (34) 

are known: A — 1, A\ = 0, and A 2 = — Ai/3. f(k + ) has support over the range — oo < 
k + < A, where A is a universal parameter in HQET which sets the deviation of 6-quark mass 
from i?-meson mass in the infinite mass limit. The values of A and Ai are not calculable 
analytically and have to be extracted from experiments or calculated on the lattice. A recent 
analysis gives (A, Ai) = (0.39 ±0.11 GeV, -0.19 q= 0.10 GeV 2 ). 

The effects of the Fermi motion of the heavy quark can be included by convoluting the 
above structure function with the differential rate [p|[7|,|T7|l , 

%= t dk + f(k + )^(mt), (35) 
aE J2E-m b dhj 

where E is the charged lepton energy in the case of semi-leptonic b —>■ u decays or the photon 

energy in the case of b — > sj decays, and dT p /dE is the rate from one loop calculation or 

the resummed version, Eq. (|20"D or Eq. (0), written as a function of the "effective mass" 

ml = rrib + k + , i.e., x = 2E/ml = MbXbI^I- The new differential rate is now a function of 

xb = IE jM.B- The addition of the structure function resums the leading-twist corrections 

and moves the endpoint of the spectrum from E = to the physical endpoint E = 

Mb/2.^ The cut rate may then be written as 



2 The true end point takes into account the final state masses, but these corrections are higher 
order in the heavy quark mass expansion. 
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2E cu t—mh 



dk + f{k + )T, 



IE, 



cut 



(36) 



where T p [2E cut / (m b + k + )] is the partonic rate with a cut at x p = 2E cut /{m b + k + ). 

VI. EXTRACTION OF \V ub \ 2 /\V t * s V tb \ 2 

From the above discussion, it is clear that a model independent extraction of IK^I 2 , 
based solely on inclusive semi-leptonic decays, will rely on a precise determinations of A and 
Ai, as well as the shape of the structure function f(k+). Knowing the first few moments of 
f{k+), A and Ai alone is perhaps not sufficient, since the shape of the structure function is 
sensitive to higher moments as well. 

Here we propose a method to extract the ratio |14b| 2 /|14*Vt&| 2 from the charged lepton 
spectrum of B — > X u £u and the photon spectrum of B — > X s 7, without invoking any 
knowledge of f(k + ). This method is based on the observation that the only difference in 
the resummations between these two processes is the g s g piece in the Sudakov exponent and 
the x v integral in the semi-leptonic decays. Applying Eq. ( |A2j ) on Eq. dl8|) , and using the 
notation introduced in Eq. (|A"l]) , we get 



M 



dx 



N 



T J dx v 6(1 - x v )x v H{x v ) (S n+1 ct n+ iJ]j +1 ^ e 



r / dx v 6(l - x v )x u H(x u ) 



Hi \ Tl 



M 



,r 

x — -N 

dx 



(37) 



where 



X ' = a s f3 \og(N + 1) « a s (3 \og(N). 



(38) 



Then taking the inverse Mellin transform of both sides and using the Convolution Theorem 
Eq. (0O|) on the right hand side we get 



dT 



dx if t 



f f 

J dx u 6x u (l — x v )H{xy) J 



m b/™„ du , dT^i 



u 



M- 1 



(39) 



u du 

where we have replaced x' with x- Let us make a few comments on these results. First, 
dT/dx and dT 1 /du are both experimentally measurable quantities. We have managed to 
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relate them in terms of theoretically known functions. The only dependence on the non- 
perturbative quantities, A and Ai, resides in the upper limit, M B jm^ for the u integral. 
Furthermore, the dependence on the structure function has been eliminated. A factor of u 2 
has also been added so that we may correctly reproduce the total rate to order of K/m b - 
Powers of u are subleading in 1/A" and thus we have the freedom to include this correction. 
However, we make no claims about reproducing corrections of order a s K/m b or (1 — a;)A/m&. 

Using the results in the Appendix the inverse Mellin transform in Eq. fl3~9|) can be easily 
performed 



M' 1 



U 



N 



9(1 -i/!i)e Js ' (aAI '/» Al . 



(40) 



Here l x / u = — log[— \og(x/u)]. Therefore 
d ^ = ^[dx 6x (l-x)*& 



M B /m b /7F7 



\v ub \ 2 
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M B /m b 



du W 



dr^ 



K 



du 



x] —z- log(l - a s p l x/u ) 
67TfJ 



(41) 



\V t * s V tb \ 2 3aC 7 (m b y J x " " du 

where the function K(x; y) in the convolution is the integral over the neutrino energy x u , 
taking into account the x dependence up to 0((1 — x) 2 ), and is given by 



K(x;y) = 6 



4a 



1 + ^(1- ^ (1) (4 + y)) 



(y + 2)(y + 3) 



3tt 
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(y + 2)2 ( y + 3f 



4a, 
3vr 



(y + 2)3 (y + 3) 3 



-3(1 -xY. 



(42) 



In the above expression we have restored the x dependence in the x v integral up to 0((1— x) 2 ) 
and ^p^(z) is the first derivative of the dig amma function 

1 d 



r(z) dz 



Viz). 



(43) 



Using the perturbative one loop level prediction for the radiative decay rate in Eq. ([II]) and 
expanding in a s leads to the one loop prediction for the semi-leptonic rate up to corrections 
of order C(a s (l - x)) and C((l - a;) 3 ). 
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The Wilson coefficient CV(/x) has been computed to next-to-leading order accuracy in 
recent years. Here we quote from Ref. |TR|T3 its value at the scale m b 



C 7 (m b ) = -0.31 + 0.48-^ + 0.03 — . (44) 

47r a. 



More explicitly, | V u ft| 2 /| V^V^,] 2 may be extracted from the relation 



\V ub \ 2 3aC 7 (m b ) 2 fi dT 

1 dx B - 



\V t * s V tb \ 2 tt(1 + 3 A/ Mb) Jx% ° dx B 



i f j r , 2 tti ' r , 

x < / dxB duBU B — — K 

Jx c n Jib CtUR 



4 

^—5- log(l - a s p Ixb/u b ) 



(45) 



ix B au B 

where we have rewritten our result in terms of the hadronically scaled variable xb = 2E/Mb, 
A = Mb — m b and the experimental cut x c B . We may use the central value for A, 390 MeV, 
with a 110 MeV standard deviation [jlfifl .P| 

One might be concerned about the existence of the Landau pole in g s ^{x)i which renders 
its inverse Mellin transform imaginary before x B — 1- In practice, this will not be a problem. 
The Landau pole in g s £ is located at (xB/uB) ma x = 1 — exp[— l/(a s Po)] ~ 0.99913, so 
numerically, one can just integrate ub in Eq. Q4"5D from x^/0.99 to 1 to avoid the Landau 
pole without incurring substantial error, since the measured rate is a smooth function. 

VII. THEORETICAL UNCERTAINTIES AND CONCLUSIONS 

By relating two disparate rates near their respective endpoints we have eliminated the 
uncertainty due to our ignorance of the structure function. However, we are left with some 
unknown uncertainty due to the assumption of quark-hadron duality. As usual, it is very 
difficult to quantify the errors due to this assumption. These errors exist in all predictions 
of heavy quark decay rates though qualitatively, we expect them to be larger in some cases 
than in others. Shifman and collaborators |2U| have argued the errors due to rotating from 



Euclidean space (this may be equated with the error in coming close to the cut in the 



3 For consistency one should use the one loop extracted value of A. 
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dispersion relation), where the OPE is performed, back to Minkowski space leads to power 
suppressed corrections. They showed this to be true in some models of QCD. We would 
expect the size of this power correction to depend upon the size of the region over which 
the differential rate is smeared. Heuristically, this just means that the partonic rate should 
approach the average of the hadronic rate as we smear over resonances. In our case, we 
smear over a region of invariant mass between M n and Md, so we would expect smearing 
over a significant number of final states, which in some sense is a measure of the sensitivity 
to parton-hadron duality. However, the spectrum is weighted towards smaller invariant mass 
states. Thus, if the decay in the end point region only went to the one pion state, we would 
expect significant errors due to duality. As the number of final states sampled increases, we 
would expect the duality errors to dwindle. 

In addition to the duality errors, there will also be the usual errors due to higher order 
corrections in l/m&. By including the factor of u 2 in Eq. fl45|) we correctly reproduce the 
total rate up to corrections which are down by A/m&. Thus the leading non-perturbative 
corrections, which have been left out, are of order A 2 /m 2 and (1— x)A/mf 1 . These corrections 
include the finite pion mass effects which shift the endpoint from Mb/2 to (M B — M 2 ) / (2Mb) ■ 
The error due to the uncertainty in A, 5A, which is dominated by statistical errors, is 
3 5A/rrib — 5% and will be reduced in the near future. 

Finally, there are the errors due to subleading terms in the exponent of the resummed 
result Eq. ([14]). This issue may be addressed by determining the effects of resummation 
which we have included. In [1(J it was found that to get a reliable prediction for B — > A s 7, 
given the present experimental cut of x c = 0.81, it was not necessary to sum the endpoint 
logs. Is the same true in the present case? We may estimate the effects of resummation 
by using a model to make a prediction for the dT 1 /dus in Eq. (^5|), and then compare 
the result for the cut rate with the rate we would find had we expanded the exp(g s i) in 
powers of a s . Doing so we find that the effects of resummation are of order ten percent. Of 
course, as emphasized in [Tt| , the true effect of resummation will depend upon the measured 



shape of the radiative end point spectrum. However, we do not expect that the effect of the 
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resummation to be so important that the inclusion of g% would be relevant. 

Before closing, we should point out that with experimental progress in reconstructing the 
neutrino momentum, it may be feasible to make a cut on the hadronic invariant mass instead 



of the lepton energy to reject charmed final states f21|| . This technique, while experimentally 
more difficult, has the advantage that the spectrum is not preferentially weighted toward 
low invariant mass states. However, as in the case discussed in this paper, the rate will again 
be sensitive to the Fermi motion of the 6-quark |22]-[24|]. Thus, an approach analogous to 



the one developed here will be necessary to get a reliable extraction of \V u b\. This work will 
be reported in a separate publication. 
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APPENDIX: 



Here we derive the inverse Mellin transform Eq. (|32|) • In Ref . |13[ the inverse Mellin trans- 
form was calculated to leading logarithm accuracy. For our purposes, we need to evaluate 
the \og(N) gi[a s log(iV)] piece to next-to-leading logarithm accuracy, whereas g2[ot s log(iV)] 
itself is already of this order. The Mellin transform pair in our case is defined as 



M[f(x);N}=r(N)= I x N - l f{x)dx 



o 



1 rC+ioo 

M- l \r(N)-x]=f{x) = — x- N f*(N)dN, (Al) 



2m JC-ioo 

We will also make use of the following two formulae which can be easily derived 

M[xf(x);N] =r(N+l), (A2) 
M- l [r(N)g*(N);x) = f f (-) g{u)— . (A3) 

Jx \Uj U 

To evaluate the inverse Mellin transform, let's begin by considering the following integral 
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C+ioo rlfJ pC'+iod A v 

S-^e^Wfi^^Wl - / ^ e *+(log*+i)Fi[a s (log*+i)] 



C-ioo iV JC'-ioo Z 



e lF 1 (cx s l) x I ^_ e z+G(lo S z,a s ,l) ^ / A4 ^ 
' Z 



where z = —N log x and 



G(logz, a s ,/) = (log 2 + I) Fi[a s {\ogz + I)) -I Fi[a s 
= \ogz(F l [a s l} + a s lF[[a s l}) 

+ log 2 z (a s F[[a s 1} + ^a s {a s l)F'([a s I] 
+ ■••, (A5) 
/= — log(— logx) , (A6) 
C = -Clogx. (A7) 

In writing down the above we have assumed that a s I is 0(1) and thus Taylor-expand with 
respect to ot s log z. The coefficient of log z in Eq. ( |A5| ) represents contributions to next-to- 
leading logarithm of the form a™l n . Terms with log n z, n > 2 only contribute to next-to- 
next-to-leading logarithm as can be seen by expanding its exponential and integrating term 
by term in z. We get precisely an expression of the form 1 + J-(a, s , I) as in Eq. (|2lf) . In 
general, if we want to perform the inverse Mellin transform to the order of a™ +k l n , we need 
to keep G up to terms of log fc+1 z and evaluate the contour integral. However, we would 
have to include in Eq. fllTf ) higher order contributions, such as a s gs, for consistency. 
Truncating G at log z order and evaluating the integral we then have 

1 j C+lo ° dz eZ zFl ( as i) +a j F l( as i) = J_ f 00 dy e C'+iy sqI + iy \Fi(a s l)+ as lF{(a s l)-l 

2ni Jc-ioo z 2n J-oo 

1 

~ T[l-F 1 (a s l)-aJF{(a s l)} ' ( 8) 

where we have changed variable from z to z = C' + iy and have used 

1 1 r°° 

= du e a+lu (a + %uY z . (A9) 

1 \z) 2ir J-oo 

Including the F 2 piece we have 
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1 r c+io ° dN 



N 



2iri Jc 



X -N e k>g(N) Fi[a a log(N)]+F 2 [ aa log(JV)] 



6(1 -x) 



J F 1 (a 3 l)+F 2 (a s l) 



where 



T[l -F 1 (a s l)-a s lFi(a s l)) 



oo fe— 1 

, k 



k=l j=0 

represents next-to-next-to-leading log contributions. 



x[l + F(a s ,l)}, (A10) 



(All) 



By taking a derivative of both sides of Eq. ( |A10| ) with respect to logx, we arrive at 

1 rC-\-ioo 

_ I dN e -N\og{x)+\og(N)F 1 [ aa \og(.N)\+F 2 [ as \og{N)\ 

2ni Jc-ioo 

x[l + F{a s ,l)}\, (A12) 



d r e l F 1 (a s l)+F 2 {a s l) 

~ X dx ~ ^ T [1 - FjaJ) - a s l F[(a s l)} 



which is relevant in processes like B — > X s 7 and Drell-Yan production. 

In the inclusive semi-leptonic decays we need to evaluate, instead, the following integral 

1 /"C+ioo ^ 



2ni Jc 



N 



. x -N e log(AT+l) Fi[a s log(7V+l)] 



dN' 



, X 1-N' e lof>(N')Fi[at a log(AT')] 



2m J N'-l 
^ 2m J N' \N' ' 



k=0 



:x) — ■ / — z e 

Z7U J Z 



Apply Eq. (|A9[) and use the following sum 
^ (-L) k 

^ r(i - 

k=0 L v 1 



F-L 



L)*e 



r(-F, -l) 



--F + k) 

where T(a,z) is the incomplete Gamma function defined as 

t a ~ l e~* dt, 

we then obtain 

1 r c+io ° dN 



r(-F) 



(A13) 



(A14) 



(A15) 



2m' Jc 



N 



_ x -N e log(JV+l)Fi[a 3 log(JV+l)] 



9(l~x)xe lFl{asl) 



£r[i-*i< 



log a;) 



a a Z) - {a s l)F[(a s l) + fc] 



0(1 -x)e 



IF^OLsl) 



logx) 



F 1 +a a lF[ 



T(-F 1 -q 8 lF{, -logx) 



(A16) 



where we have omitted the next-to- leading log contributions coming from F2, which would 
give an overall factor of exp[F 2 (a s l)]. 

On a practical note we find that, to the accuracy we are working, it is simpler to use 
large N approximation on the exponent of Eq. ( |A13| ), i.e., 

1 r c+i °° dN elog(N+1) Fl [as io g (v+i)] w J_ f c+io ° dN_ x _ N ek>s(N) Fl [as log(JV)] 

e lFi(a s l) 

0{1 ~ x) T[l-F 1 (a a l)-a a lF{(a a l)Y (A17) 



2-Ki Jc-ioo N 2ni Jc-ioo N 
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